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We propose models of two dimensional paramagnetic semiconductors where the intrinsic spin Hall 
eflFect is exactly quantized in integer units of a topological charge. The model describes a topological 
insulator in the bulk, and a "holographic metal" at the edge, where the number of extended edge 
states crossing the Fermi level is dictated by (exactly equal to) the bulk topological charge. We also 
demonstrate the spin Hall effect explicitly in terms of the spin accumulation caused by the adiabatic 
flux insertion. 
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I. INTRODUCTION 

The intrinsic spin Hall effect is a novel phenomenon in 
condensed matter physics, where a dissipationless spin 
current is proposed to be induced by an external electric 
field. The effect has been theoretically predicted both 
in p doped semiconductors with Luttinger type of spin- 
orbit (SO) coupling^ and in n doped semiconductors with 
Rashba type of SO coupling^ . After these initial propos- 
als, the issue of the stability of the intrinsic spin Hall 
effect has been intensely debated theoretically, ft is now 
broadly believed that the vertex corrections due to im- 
purity scattering exactly cancels the intrinsic spin Hall 
effect in the n type Rashba modeli^, while the vertex 
corrections vanishes for the p type Luttinger-S. and Rashba 
models^. In the latter case, the impurity scattering does 
not affect the intrinsic spin Hall effect in the clean limit. 
This conclusion is also supported by extensive numerical 
calculations^. Remarkably, the spin Hall effect has been 
experimentally observedSiiS. The experiment of Ref.^" 
was carried out in a two dimensional hole gas (2DHG) 
in the clean limit, and the effect is likely of the intrinsic 
nature. 

The next logical question in the study of the emerging 
field of the spin Hall effect concerns the dissipationless 
nature of the transport and possible quantization of the 
spin Hall conductivity. Murakami, Nagaosa and Zhang 
proposed that the intrinsic spin Hall effect can even ex- 
ist in insulators where the Fermi level lies within a band 
gap^^. In a spin Hall insulator, there is no charge current 
but spin currents, and the transport can be completely 
dissipationless. Bernevig and Zhang have proposed that 
the spin Hall effect can be quantized in two dimensioned^. 
In their proposal, the Landau levels arise from the gra- 
dient of the strain, rather than the magnetic field. 

In the present paper, we propose a new realization of 
the quantum spin Hall effect (QSHE) by specializing the 
spin Hall insulator model of Refpli to two dimensions. 
In the presence of mirror symmetry with respect to the 
xy plane, the system is shown to be a topological insu- 
lator characterized by a momentum space winding num- 



ber n g Z, with spin Hall transport carried by gapless 
edge states in a cylindrical geometry, in a way similar to 
the quantum (charge) Hall system. The evolution of the 
edge states owing to the adiabatic flux insertion can be 
traced by following Laughlin-Halperin argument for the 
integer quantum Hall effect (IQHE)iii^ii^, and it can be 
related explicitly to the spin accumulation at the bound- 
ary. Our model therefore describes a bulk topological 
insulator, and a "holographic metal" at the boundary, 
where the edge transport properties precisely encode the 
bulk topological invariant. 

The rest of this paper is organized as follows. In Sec. 
II we introduce a systematic description of the quantum 
anomalous Hall effect (QAHE) in the most general two- 
band model in two dimensions that realizes the charge 
QHE without an external magnetic field. It also pro- 
vides a helpful mathematical preparation for understand- 
ing the quantum spin Hall effect (QSHE). In Sec. Ill we 
show how the QSHE emerges in a 2-dimensional spin Hall 
insulator with an "inverted" band structure. Finally, Sec. 
IV is devoted to conclusions and discussions. 



II. QUANTIZED ANOMALOUS HALL EFFECT 

To understand the topological quantization of the spin 
Hall effect, we shall first introduce a general class of 2d 
models, called quantum anomalous Hall insulators, in 
which the charge Hall effect is topologically quantized 
in the absence of an external magnetic field. 

Historically, the first example of the QAHE was in- 
troduced by Haldane-'^^, which is a tight-binding model 
defined on a honeycomb lattice with next nearest neigh- 
bor hopping and stagger fiux. (Recently, Kane and Mele 
generalized Haldane's model and discussed the QSHE^^.) 
Similar to the usual quantum Hall effect, the QAHE is 
also a consequence of momentum space topology d'*, and 
is robust against local perturbations. We will show the 
topological nature of the spin Hall conductance in our 
general two-band model explicitly by a Kubo formula 
calculation. 
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The most general two-band Hamiltonian describing a 
2d non-interacting system can be expressed in the foUow- 
ing form: 



H = Y,H{k), H{k) = e(k) + Fd„(kK, (1) 



where tT"(a = 1,2,3) are the three Pauh matrices and 
k = {kx, ky) stands for the Bloch wavevector of the elec- 
tron. The two bands may stand for different physical 
degrees of freedom depending on the context. If they are 
the components of a spin- 1/2 electron, da(k) describe 
the spin-orbit coupling. If they correspond to the orbital 
degrees of freedoms, then ^^(k) describe the hybridiza- 
tion between bands. The discussion below is completely 
independent of the physical interpretation of the Hamil- 
tonian and leads to a general understanding for the 
conditions for the QAHE. 

The Hamiltonian |^ can be easily diagonalized to ob- 
tain the two-band energy spectrum as i?±(k) = e(k) ± 
VdCk), in which (i(k) is the norm of the 3- vector da(k). 
The Hall conductivity can be calculated using the stan- 
dard Kubo formula to be 



axy = lim —Qxy{uj + iS) , 



J- ^ tr ( J^(k)G(k, i(w„ + I/™)) 



•Jj^(k)G(k,iw„)) , 
with the current operator (i, j = x, y) 



(2) 



dki 



dki 



dki 



(3) 



and G(k, iujn) the Matsubara Green function. 

From eqs. (j2l and the Hall conductivity can be 
calculated straightforwardly. The detail of this calcula- 
tion is given in Appendix A, with the resulting Oxy given 
by 

k ^ 

where ^^(k) — (k)/(i(k) is the unit vector along 
the direction of daik)- daik) is singular if d{k) = 
^ da (k)(i"(k) vanishes for some k. However, here and 
below we are always interested in the insulating models, 
in which a full gap opens between the two bands i?+ (k) 
and £;_(k), thus £;+(k) - i;_(k) = 2Fd(k) > for aU k. 
The gap opening condition is written explicitly as 



min i?+(k) > max£^_(k) 

keBZ keBZ 



(5) 



In this case, the system becomes a bulk insulator when 
the chemical potential lies inside the gap, which implies 




FIG. 1: The Skyrmion configuration of d(k) in the Brilloune 
zone of the system JTJ with c=l,es=3.7, n = l. The vector 
d(k) starts from the north pole at the center of Brilloun zone 
and end at the south pole at the zone boundary after covering 
the unit sphere once. 



71- (k) = 1, ?T.+ (k) = for all k at zero temperature. Un- 
der such condition and taking the thermodynamic limit, 
the Hall conductivity can be simplified to 



>xy 



1 

'8^ 



dkxdkyd ■ dxd x 9yd (6) 



FBZ 



which is a topological invariant defined on the first Bril- 
lioun zone (FBZ), independent of the details of the band 
structure parametersiiS Considering d(k) : T'^ ^ as 
a mapping from the Brillioun zone to the unit sphere, 
the integrand d • 9xd x 9yd is simply the Jacobian of 
this mapping. Thus the integration over it gives the to- 
tal area of the image of the Brillioun zone on 5'^, which 
is a topological winding number with quantized value 
iirn, n e Z. Thus the conductivity axy is always quan- 
tized as axy = —n/2TT when the mapping covers n 
times. A schematic picture of a typical d(k) configura- 
tion with winding number n = 1 is shown in Fig. ^ 
Although the single-electron states in this system are 
very different from the Landau levels in the usual in- 
teger quantum Hall effect (IQHE), the quantization of 
conductivity in these two systems share the same topo- 
logical origin, which can be understood as Berry's phase 
in k-space. The exact formula plays a key role for 
the QAHE which is similar to that of the TKNN formula 
in the Landau level problem^ ^^i^? . Consequently, both of 
them are robust against weak disorder due to the topo- 
logical reason^"'^, which is well known for the IQHE case. 
The general relationship between the momentum space 
topology and the quantization of physical responses has 
been discussed extensively by Volovik in Ref.^^. 

For an explicit discussion on the QAHE and the char- 
acteristics of edge states, we consider the following choice 
of da(k) as an example: 



sin k„ 



■ sin kx , 
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dz = c {2 — cos kx — cos ky — Bs) (7) 

When V/t in Hamiltonian is large enough, the insu- 
lator condition is satisfied and the Hall conductivity 
can be shown to be 

( l/27r, < < 2 
aa^y = <^ -l/27r, 2 < e, < 4 (8) 
[ 0, > 4 or < 

where the parameter c is taken to be positive. Physically, 
this model can be understood as a tight-binding model 
describing some magnetic semiconductor with Rashba 
type SO coupling, spin dependent effective mass and a 
uniform magnetization on z direction. The experimen- 
tal realization of such QAHE will be discussed in future 
works. 

To show the behavior of edge states, one can define 
such a Hamiltonian on a strip with the periodic boundary 
condition in y-direction and open boundary condition in 
x-direction, with vanishing wave function at a; = 0, L + 1. 
In this case the Bloch wavevector ky is a good quan- 
tum number and the single-particle energy is a function 
Em{ky), (m = 1, • ■ ■ : 2L). A typical energy spectrum is 
shown in Fig. |21 For a given ky, there are 2L states, two 
of these are localized, while the rest are extended. When 
the Fermi level lies in the bulk energy gap (as represented 
in Fig. 121 by a horizontal dot line), the only gapless exci- 
tations are edge states. Similar to the usual IQHE case, 
the edge states have a definite chirality. In the present 
case, the state localized at the left edge moves with ve- 
locity Vy < and that at the right edge with Vy > 0. 
This can be seen directly from the dispersion relation of 
edge states. More generally, when axy = n/2TT, there are 
\n\ chiral edge states on each edge. 

III. QUANTUM SPIN HALL EFFECT 

Besides the interests in its own, the above discussions 
on the QAHE also serve as a natural introduction to the 
QSHE. In fact, as we shall see, the QSHE can be un- 
derstood as two copies of the QAHE; each breaks time 
reversal symmetry while the whole system remains in- 
variant under time reversal. In particular, let us consider 
the light-hole and heavy-hole bands in a semiconductor 
which can be described by the Luttinger model2^: 

^ 21 + 172 + li (s . ^ (9) 
2m m 

where the three components of S stand for the 4x4 angu- 
lar momentum matrices of a J = 3/2 electron. According 
to Ref.^^, such a 4 X 4 Hamiltonian can be reexpressed by 
introducing five Dirac F matrices F"^ — {S'*, S'^ } , a = 
1, 2, ..,5 as: 

i/(k) = e(k) + Fd„(k)F" , (10) 




-0.5 0.5 

k (in unit of 2ji) 




5 10 15 20 

lattice site coordinate i 



FIG. 2: (a) Energy spectrum of the QAHE system ^ when 
the open boundary condition is imposed in a;-direction. The 
parameters are c = l,t/V — 1/3, es — 0.5. The solid and 
dashed lines between two energy bands stand for the edge 
states at the right and left edge, respectively. The horizontal 
dotted line shows a typical chemical potential inside the gap. 
The two arrows mark the two gapless edge excitation with 
momentum ky = ±0.167r. (b) the density distribution of the 
two edge states at fermi surface, calculated for a 50x50 lattice. 

with e(k) = -^k^, V^—, 
2m m 

di — —Vikykz, d2 — —VSkxkz, 

da = -VSkxky, di = — — (k^ ~ ky) , 

and the matrices F° form the SO {5} Clifford algebra 
|pa pfcj ^ 25"''. (More details about the 50(5) repre- 
sentation of Luttinger model can be found in Appendix A 
of Ref.^'*.) By using the Clifford algebra, the Hamiltonian 
(|10|l can be diagonized to obtain the doubly degenerate 
eigenvalues: 

E±{k) = e(k) ± V^/d^dHk) ^ ^'^ ^^" k^ . (11) 

2m 

In the present work, we will focus on the spin Hall in- 
sulators described by 72 > 71 /2,^^ and specialize it to 
two dimensions. When the semiconductor described by 
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the Hamiltonian (|10|l is made into a quantum well in z- 
direction, the effective Hamiltonian can be obtained by 
adding a potential well term V{z) to the original Lut- 
tinger Hamiltonian. When Vi^z) is narrow enough, the 
system can be considered as quasi-2D in the low energy 
sector, for which one can write down a two-dimensional 
effective Hamiltonian -ff2d- The simplest way to obtain 
HiA is by replacing fc^, k\ in the original Hamiltonian H1U|) 
by their average in the lowest subband (A;^) , (fc^), respec- 
tively. When the potential is symmetric, U{z) — U{—z), 
the parity symmetry with respect to x-y plane is re- 
spected and thus (kz) = 0, which means di = 6,2 = 0, 
and the 2d Hamiltonian can be simplified to 



H2d 



e(k) + yd„(k)r",(a = 3,4,5) 
—\/3kxky, d^ = — (k^ ky) 



(12) 



2e, 



^1) 



with Cs = {k1) and k = {kx,ky). 

Noticing that Ta{a = 3,4,5) form a reducible repre- 
sentation of an 50(3) Clifford sub-algebra, it is natural 
to see the similarity of the Hamiltonian (|10l) to the two 
component one |^ proposed in the previous section. To 
see such a similarity explicitly, define 



r"''^ [^^^^]/2^, (a,5= 1,2,...,5) 



then 



[ri2, r"] = 0(a = 3, 4, 5) ^ [ri2, H2d] = O (13) 

Therefore, T^^ serves as a "conserved spin quantum num- 
ber" even in the presence of the SO coupling. The 
eigenvalues of F^^ are ±1 both with double degener- 
acy, and nad all of F" {a — 3, 4, 5) can be block- 
diagonalized simultaneously. Since they form a represen- 
tation of 5*0(3) Clifford algebra, the new expression of 
F° in the diagonal representation of F^^ can always be 
chosen as 



-/ 



(14) 



with a = 3, 4, 5, and (7^(/i = 1, 2, 3) are the Pauli matri- 
ces. In this new representation, the 2D Luttinger Hamil- 
tonian (|10|) is also block diagonal: 



Hik) 



e{-k) + Vdc,{k)a°' 



e(k) - Fd„(k)cr" 



.(15) 



In other words, the 4-component spin-3/2 system is 
equivalent to a decoupled bilayer QAHE system, each 
with Hamiltonian but with da{k.) opposite in the two 
layers. According to the definition ©, the Hall conduc- 
tivity of the two layers are opposite to each other. Since 
F^^ is odd under time reversal transformation, the two 



layers are time reversal partners, and the total Hamilto- 
nian H2d remains time-reversal invariant. The F^^-spin 
current is given by 



J = J' 



J- 



(16) 



in which jf^ is the current of electrons with F^^ = ±1. 
Thus the Hall conductivity of = a^yEy is quantized 



as 



'xy 



xy 

n 

TT 



An' 



— 2 / / dkxdkyd ■ (9xd x 9yd 



FBZ 



(n e Z) 



(17) 



In general, such a quantized Hall conductivity of the 
conserved charge F^^ leads to a non-vanishing spin Hall 
effect in the 2d insulator system Ill2fl . which is consistent 
with the three dimensional spin Hall insulator model. 
What's more, the spin Hall transport in the two dimen- 
sional system can be understood better by studying the 
edge states, as in the QHE and QAHE. To see the pic- 
ture more clearly, a tight-binding regularization of da (k) 
is specified as 



dsik) = — \/3 sin /ca: sin fcj, , 
d4{k) = VS (cosfcj, — cosfcy) , 
dzik) = 2 — Cs — coskx — cosk. 



V ' 



(18) 



which reduces to the continuum form in ean. (|12|l when 



kx,ky 0. Direct calculations show that 



2/7r, < < 4 
0, e, > 4 or e, < 



(19) 



thus the topological charge is 2 when < es < 4 and 
t/V is small. The topological charge in this system is 
larger than the previous QAHE example Q by one unit, 
since the d wave functions here "winds around" in the 
momentum space more than the p wave functions in the 
previous QAHE example. In this system there are four 
edge states on each boundary. For the F^^ = -fl(— 1) 
states, the Vy > state is localized on the left (right) 
edge, while the Vy <{) state is localized on the right (left) 
edge. The energy spectrum and the schematic diagram 
of the edge states are shown in Fig. |3| 

To study the evolution of the edge states in an infinites- 
imal electric field, we consider the Laughlin-Halperin 
gauge argumenlii^ii^. The system with the periodic 
boundary condition in y-direction and open boundary 
condition in a;-dircction can be considered as a cylinder. 
When the flux <f>(t) threading the cylinder is adiabati- 
cally turned on from $(0) = to $(r) = 27r, the electric 
field in y-direction is given by 



Ey{t) = - 



dt 



1 d^{t) 
L dt 



(20) 



For convenience, the direction of the flux <I> is chosen 
so that Ey > when $ increases. The effect of flux 
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(in unit of 2n) 



(b) 



the ground state of the system is given by |G) = 
nB„fc<^ Imfc). When a 27r flux is threaded through 
the cylinder, the final state is obtained by a trans- 
lation of Fermi sea in momentum space: \G') = 



n 



,k + 2'!T/L) Since the bulk part of \G) is a 
product of all fc's, it does not change under such a trans- 
lation, which implies that the only difference between \ G) 
and |G") occurs to the edge states near the Fermi level. 
As shown in Fig. |3l(a) by solid and hollow circles near the 
fermi level, each edge state on the Fermi surface with ve- 
locity Vy > will move out of the fermi sea and becomes 
a particle excitation since SE ~ VySk = 2TrVy/L > 0, 
while each one with Vy < will move into the fermi sea 
and leads to a hole excitation. Consequently, the final 
state |G") can be expressed as a particle-hole excitation 
state as 





=+1 




=-1 



X <: 



1 



FIG. 3: (a) The energy spectrum in the case 11811 with t/V — 4 
and Ba = 0.5. The isolated solid lines stand for the (doubly- 
degenerate) edge states and the dashed line indicates a typ- 
ical in-gap fermi energy (with /x = —4.2t). Each crossing 
of the fermi energy and the edge state spectrum defines two 
edge states on left and right boundary with opposite value of 
r'^^. The solid and hollow circles near a fermi level mark the 
particle and hole edge excitations induced by adiabatic flux 
insertion. (See the text for details.) (b)Schematic picture of 
the edge states. Each red (blue) line stands for an edge state 
with = The double arrow shows the pseudo- 

spin orientation of the current carried by the corresponding 
edge state when an electric field is applied in y-direction. For 
simplicity, only one edge state with definite F^^ eigenvalue on 
each edge is drawn. 



Ay in 



threading can be expressed by replacing ky ky 
the Hamiltonian, or equivalently, as a twisted boundary 
condition 'tplx^y + L) = e^'^il'{x,y). In such a picture, 
each single particle eigenstate (l)mky (a;, y) = u^ky (a;)e'''«^ 
will be adiabatically transformed into 4>mky{x,y){t) = 
Um,ky-Ay{x)e'^^''^~^^'^'^. In particular, when the flux 
reaches 27r, the adiabatical evolution will result in 



27r 

m,ky + — 



(21) 



\G') =\{4lAr-^^l^^^R+\G) 



(22) 



1=1 



When the Fermi level lies in the bulk energy gap. 



in which the label ± stands for the eigenvalue of F^^ 
carried by the edge state and L, R refers to the edge 
states on the left and right edge, respectively, n is the 
bulk topological number. In obtaining eqn. (|22|) . we have 
used the chirality of the edge states, i.e. vl+ > 0, wl_ < 
0,vb+ < 0,i;fl,_ > 0. (A similar analysis of the usual 
IQHE case can be found in RefiS^.) 

From Eq. (|22|l it's clear that the net effect of adiabat- 
ically turning on a 27r flux is to transfer edge states with 
F^^ = 1 from the right edge to the left one and to trans- 
fer edge states with F^-^ = — 1 in the opposite way. This 
leads to an accumulation of F^^-"spin" on the bound- 
ary. Since F^^ jg related to S"^ by S*^ = -^F^^ - F^^ as 
shown in Refi24, such an accumulation of F^^ in general 
leads to a non-vanishing spin density on the bound- 
ary. On the other hand, such an accumulation can also 
be considered as a consequence of the spin Hall current 
jx induced by the electric field Ey in Eq. H2U|) , which im- 
plies that the physically observed spin Hall conductivity is 
proportional to the amplitude of spin accumulation after 
27r-fiux threading. Since (S*^) = -^{T^"^) - (F^^), the 
corresponding spin Hall conductivity also consists of two 
parts, where the conserved part ai^^y — ^<J^y corresponds 
to a transport of F^^ spin carried by the motion of edge 
states; while the non-conserved part ax^'' is just a pre- 
cession effect due to the non-conserved nature of spin as 
represented by (F'^^) of each edge state. Consequently, it 

(c) 

is only axy that counts true transport of quantum states 
in the system and is protected by the bulk topological 
structure. These considerations give the physical justifi- 
cation of the conserved spin current operator defined in 
Ref.24. 

Finally, we consider the effect of breaking the z-axis 
mirror symmetry, which can be induced by adding an 
asymmetric potential U{z), such that U{z) ^ U{—z), in 
the Hamiltonian (|10|) . Consequently, the average c = 
(kz) in the lowest 2-d subband becomes finite. And thus 
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0' 

|i (arb. unit) 



FIG. 4: The schematic curve of conductance G versus fermi 
level /i for a ballistic quantum well in the quantum spin Hall 
regime. The plateau with Grcs = ^ shows the residual con- 
ductance contributed by 4 edge states on two edges. 

an extra term, 

= F(di(k)ri + d2(k)r2) 

= -\/3cy (fc^r^ + k^V'^) , (23) 

should be added to the two-dimensional Hamiltonian 
Since {F^^, i7a} — 0, such a term will lead to a flip 
between the states with opposite F"'^^ pseudospin. Espe- 
cially, a mixing between the left and right moving edge 
states will be induced, and thus a gap E^dgc oc cV is 
open on each edge. Consequently, the system becomes a 
fully gapped insulator when the chemical potential fi lies 
within the edge gap; however the gapless edge excitations 
still exist if /i is not in the edge gap but remains within 
the bulk gap. In this case, the spin Hall effect carried 
by the edge states can still survive, but not as robust 
as in the fully symmetric case, since it is not completely 
topology-protected . 

Spin currents in this spin Hall insulator model has also 
been discussed in Ref.^. 



IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have proposed a topological mecha- 
nism for the quantum spin Hall effect, which are realized 
in a class of two-dimensional spin Hall insulators with 
mirror symmetry. Mathematically, such QSHE is the 
cousin of the quantum Hall effect in band insulators with- 
out a magnetic field, the so-called quantum anomalous 
Hall effect (QAHE). By carefully studying the Laughlin- 
Halperin gauge argument, insights are gained for the 
physical mechanism of the quantum spin Hall transport 
carried by the edge states. It also provides an under- 
standing of the physical meaning of spin transport in the 
absence of spin conservation. 

The QSHE models discussed in this paper can be ex- 
perimentally realized in two classes of 2D semiconduc- 
tors. One class is the (distorted) zero-gap semiconductors 



such as HgTe, HgSe, /9-HgS and a-Sn. The other class 
is the narrow-gap semiconductors such as PbTe, PbSe 
and PbS— . As proposed above, topological quantiza- 
tion of the spin Hall effect shows up in the cases with 
mirror symmetry with respect to the z-axis, realizable 
when the 2D material is trapped in a symmetric quantum 
well. Once the quantum spin Hall effect is realized in the 
ground state, it is protected against thermal fluctuations 
by the bulk energy gap Eg. At temperature T <^ Eg/ks, 
the quantum spin Hall effect is expected to be observ- 
able. Take HgTe/Hgi_xCdxTe (001) quantum well for 
an example. As calculated in Ref^^, the gap Eg between 
LH and HH bands is of order lOmeV for x = 0.7, which 
means the quantum spin Hall effect can be observed in a 
wide temperature range T ^ \QQK. Once the spin Hall 
effect in the pure system is established, it is not signifi- 
cantly dependent on the mobility of the material. 

However, the magnitude of edge spin accumulation in 
a steady state is in general dependent on the spin re- 
laxation mechanism and the disorder in the system. To 
avoid ambiguity in estimating the edge spin accumula- 
tion, here we propose a more definite experimental pre- 
diction of the QSHE. As is shown above, there are two 
pairs of gapless edge states on each boundary, which are 
chiral in spin transport but nonchiral in the charge chan- 
nel. In other words, each pair of edge states (with oppo- 
site F^^ on the same boundary) is equivalent to a spinless 
Luttinger liquid in two-terminal measurements of charge 
conductance. Consequently, there are in total 4 pairs of 
edge states in the system, which will make a quantized 
contribution^^. 



to the longitudinal charge conductance when the system 
is in a ballistic regime. The schematic curve of longitu- 
dinal conductance versus the Fermi level is shown in Fig. 
^ Compared to the vanishing conductance in a trivial 
insulator, such a residual conductance provides a simple 
probe of the topologically non-trivial edge states in the 
QSHE. 
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APPENDIX A: KUBO FORMULA 
CALCULATION OF a.^y 

Here we propose a derivation of eqn. © from the Kubo 
formula. The single particle Green function correspond- 
ing to the Hamiltonian ^ can be written as 

G{k,iujn) = (iw„ - -ff (k))""^ 
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^+ I ^- (Al) Then the charge Hall conductivity can be calculated us- 

-E+{k) iujn- E_{k) ' ^ ing the Kubo formula @: 



with P± = (l±d„(k)CT' 



ixy 



{iv,n) = -i-^tr(J^(k)G(k,i(w„ + j/„))Jy(k)G(k,zw„)) 

k,n 

^ ^ V- tr(J,(k)P,(k)J^(k)Ft(k)) 

^^^^ (* K + ^m)- Esik)) {ttOn ~ Et{k)) 



r ^ r^ I . ■x^ * tr ( J,(k)P,(k) Jy(k)Ft(k)) 

hm — (3j;„(tj + id) = -— > > — ^ — — — „ [ntik) - n^fk)) 

= E ""'whiikY^^ [tr (J.(k)P+(k)J,(k)P_(k)) - fe.c] . (A2) 

k 

Here E±{k) ~ e(k) ± Vd{k) is introduced in the last line. The trace in eqn. I|A2|I can be worked out by substituting 
Ja;,y(k) and P±i)i.) with their definition and ljAl|l . respectively: 



n_(k) — rt+(k) 
(A3) 



r 



For the last step, it should be noticed that only the three- coupling V , as expected from topological considerations, 
cr" terms in the expansion make non- vanish contribu- This finishes our derivation of the formula Q for the spin 
tion to the trace. Thus the final result is independent of Hall conductance. 
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